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Genetic linkage data have shown that alterations of the BRCA I gene are responsible for the majority of hereditary breast and 
ovarian cancers. BRCAI germline mutations, however, are found less frequently than expected. Mutation detection strategies, 
which are generally based on the polymerase chain reaction, therefore focus on point and small gene alterations. These 
approaches do not allow for the detection of large gene rearrangements, which also can be involved in BRCAI alterations. 
Indeed, a few of them, spread over the entire BRCAI gene, have been detected recently by Southern blotting or transcript 
analysis. We have developed an alternative strategy allowing a panoramic view of the BRCAI gene, based on dynamic molecular 
combing and the design of a full four-color bar code of the BRCAI region. The strategy was tested with the study of four large 
BRCAI rearrangements previously reported. In addition, when screening a series of 10 breast and ovarian cancer families 
negatively tested for point mutation in BRCA 1/2, we found an unreported 17-kb BRCAI duplication encompassing exons 3 to 
8. The detection of rearrangements as small as 2 to 6 kb with respect to the normal size of the studied fragment is achieved 
when the BRCAI region is divided into 10 fragments. In addition, as the BRCAI bar code is a morphologic approach, the direct 
observation of complex and likely underreported rearrangements, such as inversions and insertions, becomes possible. 
©2001 Wiley-Uss, Inc. 



INTRODUCTION 

Genetic predisposition associated with a high 
risk of breast and ovarian cancer is one of the most 
common pathologic genetic conditions, with a car- 
rier frequency estimated at 1/500. Alterations of 
either the BRCAI or BRCA2 gene, transmitted in 
an autosomal dominant mode, account for the ma- 
jority of familial breast and ovarian cancer cases 
(Ford et al., 1998). Indeed, linkage analyses per- 
formed on a set of 94 families with at least four 
cases of breast cancer and at least one ovarian 
cancer case have shown that BRCAI and BRCA2 are 
involved in 95% of families, BRCAI being altered 
in 80% of such cases (Ford et al., 1998). Although 
numerous efforts in screening for BRCAI have 
been made, only two-thirds of expected mutations 
are detected at present (Gayther and Ponder, 1997; 
Stoppa-Lyonnet et al., 1997; Ford et al., 1998). 
Most reported approaches (single-strand conforma- 
tion polymorphism, heteroduplex analysis, dena- 
turing gradient gel electrophoresis, protein trunca- 
tion test, direct sequencing) are based on 
polymerase chain reaction (PCR) techniques and 



therefore focus on point and small mutations (Gay- 
ther and Ponder, 1997). 

The low BRCAI mutation detection rate may be 
due in part to the existence of large rearrangements 
not detected by such approaches. Indeed, PCR- 
based methods can miss an altered allele when it is 
hidden by the persisting normal allele. Supporting 
this hypothesis, 15 different large rearrangements 
have been detected recently by Southern blotting 
or analysis of BRCAI lymphocyte transcripts 
(Payne and King, 1997; Petrij-Bosch et al., 1997; 
Puget et al, 1997, 1999a, b; Swensen et al., 1997; 
Carson et al., 1999; Montagna et al., 1999; Rohlfs et 
al., 1999, 2000). The characterized rearrangements 
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are due to unequal recombination events between 
Alu sequences that cover 41.5% of BRCA1 introns 
(Smith et al., 1996). First estimations of the fre- 
quency of large rearrangements in the BRCA1 mu- 
tation spectrum range from 15% to 36%, the high- 
est value reflecting founder effects in the Dutch 
population (Petrij-Bosch et al., 1997; Puget et al., 
1999a). 

The reported large rearrangements range from 
0.5 to 23.8 kb and are spread over the entire BRCA1 
gene. In addition, three of the reported deletions 
remove portions of the 5' end of BRCA1 and its 
surrounding regions (Swensen et al., 1997; Puget et 
al., 1999a). Thus, screening BRCA1 for large rear- 
rangements requires a panoramic view of the 81 kb 
of the gene and its flanking regions. With this goal, 
we have chosen the dynamic molecular combing 
technique associated with fluorescence in situ hy- 
bridization (FISH) (Bensimon et al., 1994; Micha- 
let et al., 1997). Combing relies on homogeneous 
stretching of DNA molecules, which is achieved 
using a simple apparatus that dips a silanized cov- 
erslip into DNA solution and pulls it out at a 
constant vertical speed to obtain DNA molecules 
aligned and stretched at a constant rate of 2 kb/jxm. 
We have developed a FISH strategy inspired by 
Florijn et al. (1995) and Michalet et al. (1997). A 
complex BRCA1 full four-color bar code was there- 
fore designed, allowing for analysis of the whole 
region divided into 10 fragments. 

Here we report the usefulness of this approach 
by analyzing four BRCA1 rearrangements previ- 
ously described (Montagna et al., 1999; Puget et al., 
1999a, b). In addition, the efficiency of the ap- 
proach was examined on a series of 10 breast and 
ovarian cancer cases that tested negative for the 
BRCA1/2 point mutation. 

MATERIALS AND METHODS 

Preparation of Genomic DNA for Combing 

Total human genomic DNA was obtained from 
peripheral blood lymphocytes or lymphoblastoid 
cell lines. Because use of high-molecular-weight 
DNA is a limiting factor, cells were embedded in 
low-melting agarose blocks, as described previ- 
ously (10 6 cells per block, equivalent to 6 u,g of 
DNA) (Michalet et al., 1997). DNA was diluted to 
approximately 1.5 jxg/mL in 400 mM MES-NaOH 
at pH 5.5 (Sigma, St. Louis, MO). Combing was 
done using the Molecular Combing Apparatus™ 
(Institut Pasteur, Paris, France) and with silanized 
surfaces (22 X 22 mm) (Michalet et al., 1997). 



Labeling of Probes 

Two types of probes were used: genomic DNA 
cloned in P, -derived artificial chromosomes (PACs) 
103014 and 44B1 (from chromosome 17 library, 
LLNL, Livermore, CA) (Brown et al., 1995) or 
cosmid ICRFcl05D06121 (from RZPD, Berlin, 
Germany) (Brown et al., 1996) and products of 
long-range PCR (LR). PACs and the cosmid were 
grown and DNA extracted with the very low copy 
cosmid purification protocol (Qiagen, Valencia, 
CA). Long-range PCR products were generated 
from PACs with the Expand™ long template PCR 
system (Roche Diagnostics, Mannheim, Germany). 
Primers used are available on the Institut Curie 
Web site (http://www.curie.net/genetique) or upon 
request. Each probe (1 fig) was labeled by random 
priming (BIOPRIME DNA Labeling System; 
GIBCO, Gaithersburg, MD): PAC 103014 with fluo- 
rescein- 11 -dUTP (FluoroGreen; Amersham, Piscat- 
away, NJ), cosmid with digoxigenin-ll-dUTP (Roche 
Diagnostics), LR9-12 (6,5 kb) and LR24-3' (10 kb) 
with biotin-14-dCTP (GIBCO), and LR16-19 (8 kb) 
with both bio tin and digoxigenin. Each probe was 
purified with the QIAquick PCR purification kit 
(Qiagen). Nine hundred nanograms of each probe 
was ethanol-precipitated with 12X excess human 
Cot-1 DNA (1 mg/mL; GIBCO). 

Hybridization and Probe Detection 

Hybridization was performed as previously de- 
scribed (Michalet et al., 1997). Biotinylated, digoxi- 
genin- and fluorescein-labeled probes were de- 
tected with AMCA, Texas Red, and fluorescein, 
respectively. Four or five successive layers of anti- 
bodies were used as follows: (1) and (3) AMCA- 
avidin 1/10 (Vector Laboratories, Burlingame, CA) 
plus anti-dig-Texas Red (mouse) 1/50 (Jackson 
Immunoresearch, West Grove, PA) plus anti-fluo- 
rescein (rabbit) 1/10 (Molecular Probes, Eugene, 
OR); (2) and (4) biotinylated anti-avidin 1/50 (Vec- 
tor Laboratories) plus anti-mouse- dig 1/10 (Roche 
Diagnostics) plus anti-rabbit-fluorescein (goat) 
1/10 (Vector Laboratories); and (5) AMCA-avidin 
1/10 plus anti-dig-Texas Red (mouse) 1/50. 
Washes and antibody incubations were carried out 
as previously described (Michalet et al., 1997). 

Image Analyses 

Full signals, corresponding to the BRCA1 gene, 
were observed under an epifluorescence Leica 
DMRB microscope and captured with IPLab Spec- 
trum-SU2 software (Vysis, Downers Grove, IL) us- 
ing a NU 200 CCD camera (Photometries, Tucson, 
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AZ). Image analyses were performed with Car- 
tographiX software (X. Michalet, Institut Pasteur), 
which converts fluorescent spots into intensity 
curves. Measurements of each fragment in a series 
of signals are reported on each histogram, resulting 
in an estimation of the fragment size with its stan- 
dard deviation due to the constant stretching rate 
of 2 kb/u,m. CartographiX draws normal curves, 
thereby identifying two putative allele populations. 

RESULTS 

The BRCA I Bar Code Design 

The genomic structure of the BRCAI region and 
the probes used are shown in Figure la. Counter- 
staining of the BRCAI region was performed with a 
PAC. A cosmid was used as a 5' probe and a 
long-range PCR product, LR24-3', at the 3' end. 
Since we considered the interval between these 
probes to be too large (76.25 kb) (Smith et al., 
1996) to detect rearrangements of less than 5 kb — 
and to improve their localization — two additional 
LR products, LR9-12 and LR16-19, were hybrid- 
ized (Fig. la). By using combinations of different 
fluorescences of probes that divided the BRCAI 
region into 10 fragments, the following full four- 
color bar code was expected: [NBRl]-cosmid- 
[NBR2]-cosmid-[intron 2-intron 8]-LR9-12-[in- 
tron 12-intron 15]-LR16-19-[intron 19-exon 24]- 
LR24-3'. 

BRCAI Bar Code of Two DNA Controls 

DNA samples from IC918 and IC1063 families, 
each harboring a point mutation in BRCAI or 
BRCA2 (Stoppa-Lyonnet, unpublished data), were 
used as controls to optimize the BRCAI bar code 
strategy. The DNA samples were combed and si- 
multaneously hybridized with the five probes. Af- 
ter screening each slide, full signals of the BRCAI 
gene were compiled and organized with the same 
orientation, to facilitate qualitative and quantita- 
tive analyses. Two IC918 full signals are presented 
in Figure lb. Hybridization of the five probes was 
observed, leading to the expected bar code of the 
BRCAI region. 

The BRCAI gene was divided for measurement 
into seven fragments ranging from 6.5 to 24 kb. 
Signals of each fragment were then measured by 
CartographiX software. The seven IC918 histo- 
grams are shown in Figure lc. A single peak with a 
normal shape was seen for each histogram. The 
expected and measured sizes (averages and stan- 
dard deviations) of the seven fragments are re- 
ported in Table 1. The observed sizes of LR sig- 



nals were close to their expected sizes, with 
standard deviations ranging from 0.4 to 1.6 kb. For 
green PAG signals, the observed sizes were lower 
than the expected sizes deduced from the genomic 
sequence (Smith et al., 1996). For these larger 
fragments, standard deviations ranged from 2.0 to 
3.0 kb. The variability of a measured fragment is 
defined as the standard deviation divided by aver- 
age size. For both controls, the observed variations 
range from 6% to 18% (average, 13%; Table 1). 
Variability may be linked to the beaded aspect of 
signals likely due to the lack of hybridization signal 
on repeated sequences and to the mild heteroge- 
neity of the coverslip treatment (Michalet et al., 
1997). Because the variability is independent of 
fragment size, however, the smaller the studied 
fragment, the higher the precision of its measure- 
ment (Table 1). 

Finally, when studying a rearrangement, the two 
fragments corresponding to the two allele popula- 
tions are expected to be detected through the iden- 
tification of two peaks on a histogram, discrimi- 
nated by twice the standard deviation (2SD). 
Considering the 6.5 kb of LR9-12 and the 24 kb of 
the largest green PAC fragment, [intron 2-intron 8], 
deletions in these fragments of at least 1.7 kb and 
6.2 kb, respectively, should be visible. According to 
the experimental protocol, 100 equivalent genomes 
per slide are expected, a quantity sufficient for 
thorough analysis. Indeed, mean fragment sizes for 
control measurements are similar for 20, 40, or 70 
signals. If the quantity of signals does not appear to 
be limiting, the quality of DNA is crucial to avoid 
broken signals. 

Detection of a Large Deletion Removing One 
of the Probes 

Puget et al have identified by Southern blotting 
a 23.8-kb deletion that removes exons 8 to 13 of 
BRCAI in the IC568 family (Puget et al, 1999a). 
After combing and hybridization, this large BRCAI 
deletion was detected easily with few fields of 
view. Some full signals without LR9-12 were ob- 
served, thereby allowing us to distinguish between 
the two alleles (Fig. 2a). To evaluate the size of the 
deletion better, the region between cosmid and 
LR16-19, which includes both breakpoints, was 
measured on the entire series of signals (Fig. 2b; 
Table 1). As expected, two peaks were obtained, 
representing the two allele populations. The size 
and the standard deviation of the deletion were 
therefore estimated as the difference between the 
two averages: 20.0 ± 9.6 kb. The known size of 
23.8 kb falls within this interval. The large standard 
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deviation is linked to the size of the analyzed 
fragment (expected to be 50 kb for the normal 
allele and 26.4 kb for the deleted allele). From 
these results, it appears that when a fluorescent 
fragment is removed, the color bar code strategy 
gives qualitative information: the existence of large 
BRCAI rearrangements, their locations, and the 
estimation of their minimal sizes. Nevertheless, it 
remains difficult to evaluate rearrangement size 
precisely. 

Detection of a 6-kb Duplication 

A duplication, recently identified by lymphocyte 
transcript analysis and precisely characterized by 
PCR, was shown to be a 6-kb duplication of exon 
13 and its surrounding regions (Puget et al., 1999b). 
One of the DNA samples harboring this duplica- 
tion, IARC3653, was analyzed. Through micro- 
scope screening and by signal alignment, it was 
difficult to distinguish two allele populations (data 
not shown). Measurements of the fragment of in- 
terest, [intron 12-intron 15], however, yielded two 
peaks (Fig. 3a; Table 1). The size of the duplica- 
tion was estimated at 5.8 ±1.8 kb, compared with 
the expected 6 kb. This result confirms the reso- 
lution capacity of this approach of at least 2SD, 
which was expected to be 5.1 kb for this 19.6-kb 
fragment. 

Detection of Two 3-kb Deletions 

A 3-kb deletion, which removes exon 17, was 
identified in the Italian family UP-B74 (Montagna 
et al., 1999). When UP-B74 signals were aligned, it 
was possible to distinguish several full signals with 
a shortened LR16-19 probe (Fig. 3b). Measure- 
ments of this fragment showed two peaks puta- 
tively associated with two allele populations and 



with average lengths differentiated by 3.2 ± 0.8 kb 
(Fig. 3c; Table 1). We expected to detect the 
UP-B74 deletion, since it is located in a fragment of 
8 kb with a 2SD of 2.1 kb. This value is in agree- 
ment with our estimation of the resolution for the 
detection of rearrangements that fall within one 
probe. 

To determine whether we are able to detect a 
rearrangement of a size less than 2SD, we exam- 
ined the IC788 DNA carrying a 3-kb deletion re- 
moving exon 15 within the 19.6-kb fragment [in- 
tron 12-intron 15] (Puget et al., 1999a). 
Measurements of the fragment of interest showed 
two peaks, leading to the estimation of the size of 
the deletion at 3.6 kb; however, there was a high 
standard deviation of 2.8 kb (Fig. 3d; Table 1). 
Whereas the detection of a 3-kb deletion in an 8-kb 
fragment is feasible, that of a 3-kb deletion in a 
20-kb fragment is at the limit of the resolution of 
this approach. 

Identification of an Unreported BRCAI Duplication 
by Color Bar Code Screening 

The color bar code strategy was applied in the 
analysis of a series of 10 French breast and ovarian 
cancer families that tested negative for a point 
mutation in BRCAI 12. IC827 DNA showed several 
full signals larger than the others. Through align- 
ment, the enlarged region appeared to be located 
between cosmid and LR9-12 (Fig. 4a). Measure- 
ments of this region gave two peaks, suggesting the 
existence of at least a 10-kb duplication in one 
allele (Fig. 4b; Table 1). To confirm this observa- 
tion, cosmid and LR9-12 were substituted with two 
additional probes: LR2-3 and LR5-8 (Fig. 4a). The 
region between exons 3 and 8 appeared to be 
duplicated in tandem, as shown by the chimeric 



Figure I . a: The BRCAt region and the probes used for bar coding of the region. The BftCA I gene is spread over 81 kb (Smith et a)., 1 996) and has 
a common promoter with the NBR2 gene (next to BRCAI gene 2), which covers approximately 30 kb (Xu et al., 1997). The BRCAI pseudogene 
(VBRCAI) lies next to NBR2 and corresponds to a partial duplication of BRCAI, from the promoter region to intron 2 (Barker et al., 1996; Brown et 
al., 1996). The NBRI gene (Next to BRCAI gene /) is located 5' to VBRCAI and may be phylogenetically at the origin of NBR2 (Campbell et al., 1994; 
Xu et al., 1997). The figure is not drawn to scale and is adapted from Brown et al. (1996) and Xu et al. (1997). Probes used for the full four-color bar 
code of BRCAI: the PAC I03OI4 insert covers the region between the first exons of NBR/ and the 3' BRCAI UTR (exon 24), which is approximately 
120 kb long (Brown et al., 1995, 1996; Smith et al., 1996; Xu et al., 1997). By fluorescence in situ hybridization, this clone hybridized on chromosome 
17 and also on chromosome 4, suggesting that this PAC is chimeric. The cosmid ICRFcl05D06l2l clone used in our experiment is rearranged, since 
the insert was found to be less than 10 kb (Brown et al., 1 996; Xu et al., 1997). By sequencing, one end was localized at the beginning of BRCAI intron 
2. Considering the size of the insert, the other end was expected to lie in intron I of NBR2. Thus, our recombinant clone does not spread over the 
entire NBR2 gene, as previously described (Brown et al., 1996; Xu et al., 1997). In addition, our clone was expected to cross-hybridize with VBRCAI. 
The PAC 44BI insert has been located in the 3' BRCAI region (Brown et al., 1995). Its BRCAI end lies in intron 13, and the other end is expected to 
be located upstream of the BRCAI 3' end (Smith et al., 1996). Because this PAC appeared to be too large to be used as a 3' probe, a long-range PCR 
product, LR24-3', was generated from this clone, covering BRCA I exon 24 to 10 kb downstream. Two other probes, generated by long-range PCR from 
PAC I03OI4, were added: LR9-I2, from exons 9 to 12 covering 6.5 kb, and LRI 6-19, from exons 16 to 19 covering 8 kb. b: IC9I8 control DNA: two 
full signals showing the bar code of BRCAI (IC: Institut Curie). After screening the IC9I8 slide, 74 full signals were captured. To facilitate their view 
and analyses, full signals were aligned. All probes hybridized well. Since the four shortest probes overlap with the counterstaining PAC, color mixes 
should have been obtained at their locations. The lack of PAC hybridization signal in these overlapping regions is probably due to its low molarity 
compared with that of the four other probes, since the same quantity was used for each of the five probes, c: The seven IC9 1 8 histograms, with microns 
in abscissa and the number of signals in ordinate. The analyzed fragment and its expected size are reported in the top part of each histogram. The 
observed sizes and their standard deviations are indicated in kilobases below each histogram. 
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deleted peak 24.8 ± 4.0 kb 

normal peak 44.8 ± 5.6 kb 

■► deletion 20.0 ± 9.6 kb 

Figure 2. IC568 DNA harboring a 23.8-kb deletion (from exons 8 to nature of the rearrangement, its localization, and minimal size, b: His- 

1 3) (Puget et a!., 1 999a). a: Two normal full signals ( I ) and two full signals togram of measurements between cosmid and LR 1 6- 1 9 on both normal 

without LR9-I2 probe (2). Sixty-two full signals were obtained — 30 and deleted signals. Two peaks representing the two allele populations 

signals with LR9-I2 and 32 without — allowing the distinction between are detectable. The size of the deletion is obtained by determining the 

the two allele populations and yielding qualitative information about the difference between the two peaks. 



hybridization signal of LR5-8 and a small part of 
LR2-3. Fine characterization using long-range 
PCR in the region of interest showed a 17,247 bp 
duplication of exons 3 to 8 (nucleotides 11967- 
29213 [GenBank accession number L78833)), 
which is the result of a homologous recombination 
between an Sq Alu in intron 2 and an Sp Alu 
sequence in intron 8, which share 83% homology 
(data not shown). Under the microscope, the other 
DNA samples did not seem to have large rear- 
rangements, and signal measurements did not de- 
tect additional rearrangements (data not shown). 

DISCUSSION 

We have shown that the strategy combining 
DNA combing and full four-color bar coding is 
useful for detecting large BRCA1 rearrangements, 
since a panoramic view of the gene and its sur- 
rounding regions can be obtained. The detection of 
rearrangements as small as 2 to 6 kb with respect to 
the normal size of the studied fragment is achieved 



when the BRCA1 region is divided into 10 frag- 
ments. The efficiency of our approach was verified 
with the detection of a previously unreported 
BRCA1 duplication when screening 10 breast and 
ovarian cancer cases that had tested negative for 
the B RCA 1/2 point mutation. This bar code is ex- 
pected to be able to detect 80% of the large rear- 
rangements previously reported in the literature. 
Hybridization with additional probes and use of 
other colors, splitting the BRCA1 gene into frag- 
ments shorter than 10 kb, should lead to the de- 
tection of rearrangements smaller than 2 kb. An 
optimized bar code and the analysis of a larger 
series would allow for accurate estimation of the 
contribution of large rearrangements in the BRCA1 
mutation spectrum. 

Detection of large rearrangements represents a 
recurrent challenge in molecular pathology. South- 
ern blotting, which was the first tool used to detect 
large rearrangements, relies mainly on the detec- 
tion of abnormal restriction endonuclease frag- 
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Figure 3. a; IARC3653 DNA harboring a 6-kb duplication (exon 13; IARC: International Agency for 
Research on Cancer) (Puget et al., 1999b). Histogram of measurements between LR9-I2 and LRI6-I9, 
[intron 12-intron 15], on 54 full signals. Two peaks are detected, b, c: UP-B74 DNA harboring a 3-kb 
deletion (exon 17; UP: University of Padova) (Montagna et al., 1999). b: Two normal full signals (I) and two 
full signals with a smaller LR 1 6- 1 9 probe (2). c: Histogram of LR 1 6- 1 9 measurements on 46 full signals. Two 
peaks are observed, d: IC788 DNA harboring a 3-kb deletion (exon 15; Puget et al., 1999a). Histogram of 
[intron 12-intron 15] measurements on 46 full signals. Two peaks are detected. 
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the region between cosmid and LR9- 1 2: [intron 2-intron 8]. (2) Hybrid- 
ization of two additional long-range products covering exons 2 to 3 of 
8.5 kb (LR2-3) and exons S to 8 of 7 kb (LR5-8), in substituting for 



cosmid and LR9-I2, respectively. The lower signal contains an additional 
chimeric hybridization pattern of probes LR5-8/LR2-3, indicating the 
duplication junction, b: Histogram of [intron 2-intron 8] measurements 
on 57 full signals. Two peaks representing the two allele populations are 
detectable. 



ments induced by rearrangements. This strategy 
may fail when the target sequence recognized by 
the probe also is deleted. The difficulty encoun- 
tered in distinguishing between the presence of 
one or two copies in the case of a deletion is 
augmented when distinguishing between two or 
three copies in the case of a duplication. A second 
tool is the analysis of transcript size in the available 
tissue (lymphocytes in most cases). It requires con- 
ditions that are not necessarily fulfilled, however: 
transcription of the studied gene and stability of 
altered allele transcripts. The approach described 
in this study may avoid these disadvantages. In 
addition, because it is a morphologic approach, the 
direct observation of complex and likely underre- 
ported rearrangements, such as inversions and in- 
sertions, becomes possible. This new strategy, un- 
til now primarily used in research laboratories, can 
be implemented for the molecular diagnosis of 
other genetic diseases as well. Indeed, the bar code 
of other genes may be designed easily using artifi- 
cial clones and/or genomic sequence information. 



Development of software allowing, on the one 
hand, the automatic capture of signals with the 
help of motorized microscope stages and, on the 
other hand, their automatic measurements would 
favor the use of this new strategy by clinical genet- 
ics laboratories. 
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